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ABSTRACT 

Neon abundances are derived in four Galactic WC stars - j 2 Vel (WR11, 
WC8+07.5III), HD 156385 (WR90, WC7), HD 192103 (WR135, WC8), and WR146 
(WC5+08) - using mid-infrared fine structure lines obtained with ISO/SWS. Stel- 
lar parameters for each star are derived using a non-LTE model atmospheric code 
OO i (Hillier & Miller 1998) together with ultraviolet (IUE), optical (INT, AAT) and in- 

frared (UKIRT, ISO) spectroscopy. In the case of 7 2 Vel, we adopt results from De 
Marco et al. (2000), who followed an identical approach. 

ISO/SWS datasets reveal the [Nem] 15.5/mi line in each of our targets, while 
[Nell] 12.8/um, [Siv] 10.5/im and [S in] 18.7/im are observed solely in 7 2 Vel. Using a 
| method updated from Barlow et al. (1988) to account for clumped winds, we derive 

"-■»« . Ne/He=3-4x 10 -3 by number, plus S/He=6x 10 -5 for j 2 Vel. Neon is highly enriched, 

*-£h ■ such that Ne/S in j 2 Vel is eight times higher than cosmic values. However, observed 

Ne/He ratios are a factor of two times lower than predictions of current evolutionary 
Q ■ models of massive stars. An imprecise mass-loss and distance were responsible for the 

| much greater discrepancy in neon content identified by Barlow et al. 

c/3 . Our sample of WC5-8 stars span a narrow range in T* (=55-71kK), with no trend 

towards higher temperature at earlier spectral type, supporting earlier results for a 
larger sample by Koesterke & Hamann (1995). Stellar luminosities range from 100,000 
to 500,0001,0, while 10~ 51 < M/{M Q yr- r ) < lO" 4 5 , adopting clumped winds, in 
\ which volume filling factors are 10%. In all cases, wind performance numbers are less 

^ • than 10, significantly lower than recent estimates. Carbon abundances span 0.08 < 

C/He < 0.25 by number, while oxygen abundances remain poorly constrained. 

Key words: stars: Wolf-Rayet - infrared - abundances - evolution - individual: 
WR11, WR90, WR135, WR146 
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1 INTRODUCTION 

Understanding the physics of massive (Mi n i t <;25 Mq) stars, 
their atmospheres, radiation, and evolution is important for 
many aspects of astrophysics since their powerful winds af- 
fect the energy and momentum balance of the interstellar 
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medium (ISM). In particular, Wolf-Rayet (WR) stars pro- 
vide crucial tests of nuclear reaction chains. 

However, quantitative analysis of such stars, repre- 
sents a formidable challenge. The assumptions of plane- 
parallel geometry and local thermodynamic equilibrium 
(LTE), which are often adopted for lower luminosity stars, 
are inadequate. Nevertheless, the properties of a large sam- 
ple of carbon sequence (WC-type) WR stars have now 
been quantitatively derived using detailed models, account- 
ing for non-LTE effects, spherical geometry and an ex- 
panding atmosphere (Howarth & Schmutz 1992; Koesterke 
& Hamann 1995). However, recent studies have demon- 
strated that clumping (Moffat et al. 1988; Hillier 1991, 1996; 
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Schmutz 1997) and line blanketing (Schmutz 1997; Hillier & 
Miller 1998) may have a significant effect on the derived 
physical properties of WR stars. 

Overall, evolutionary predictions for massive stars (e.g. 
Meynet et al. 1994) are in good agreement with the observed 
properties of Wolf-Rayet stars, van der Hucht & Olnon 
(1985) derived a Ne/He ratio for 7 2 Vel (WR11, HD 68273, 
WC8+07.5) from IRAS space-based observations, which 
was found to be in good agreement with expectations. How- 
ever, Barlow et al. (1988, hereafter BRA88) identified a nu- 
merical flaw in these calculations and added new ground- 
based observations to reveal Ne/He=1.0±0.35x 10 -3 , a fac- 
tor of six lower than predicted. Is this discrepancy due to 
failure of evolutionary models, peculiarities of the 7 2 Vel bi- 
nary system, or incorrect assumptions for the stellar prop- 
erties of the WC8 star? The combination of Short Wave- 
length Spectrometer (SWS), Infrared Space Observatory 
(ISO) observations of a larger sample of WC stars, plus re- 
cent progress in quantitative modelling of WR stars, should 
provide a definitive answer to this question. 

Recent space-based spectroscopy of WR11, obtained 
with the Infrared Space Observatory (ISO), was presented 
by van der Hucht et al. (1996). who quoted excellent agree- 
ment of the fine-structure neon line fluxes with observa- 
tions used by BRA88. Morris et al. (1998) subsequently re- 
estimated Ne/He using contemporary information on the 
stellar distance (van der Hucht et al. 1997; Schaerer et al. 
1997) and mass-loss rate (Stevens et al. 1996) of WR11, 
which revealed a considerable neon enrichment. Willis et 
al. (1997, hereafter Paper I) also used ISO/SWS to ob- 
serve WR146, another WC+O binary, again revealing sig- 
nificantly enriched neon. Recently, Morris et al. (2000) anal- 
ysed ISO/SWS observations of WR147 (WN8+OB), reveal- 
ing neon, sulphur and calcium abundances in good agree- 
ment with cosmic values, as expected by evolutionary mod- 
els. 

In this paper, we supplement results from Paper I with 
other WC stars for which ISO-SWS observations are avail- 
able, HD 156385 (WR90, WC7) and HD 192103 (WR135, 
WC8). We also re-analyse WR146 (WC5+08) using more 
sophisticated analysis techniques and in the light of other 
observational evidence (Niemela et al. 1998; Dougherty et 
al. 2000). A neon abundance is also re-derived for WR11, 
following recent spectroscopic results from De Marco et al. 
(2000). 

The outline of the present work is as follows. New UV, 
optical and IR observations of our programme stars are pre- 
sented in Sect. ^, with basic properties discussed in Sect. || 
The spectroscopic technique is introduced in Sect. ^, and 
applied in Sect. |E| Our spectroscopic results are discussed 
in Sect. ^| with neon and sulphur abundances derived in 
Sect. m. Finally, conclusions are reached in Sect. H. 



2 OBSERVATIONS 

The programme Galactic WC stars are listed in Table jj] 
where we give the various catalogue names and our adopted 
spectral types, following Smith et al. (1990) and Crowther et 
al. (1998). We will refer to our programme stars by their WR 
catalogue number (van der Hucht et al. 1981). Wind veloci- 
ties are taken from UV resonance line measurements (Prinja 



Table 1. Programme Galactic WC stars. Spectral types are ob- 
tained from Crowther et al (1998), which follows Smith et al. 
(1990) except that revised WCE criteria are adopted. Wind ve- 
locities arc taken from (a) Prinja et al. (1990), (b) Willis et al. 
(1997), or (c) St Louis et al. (1993) 



WR 


HD 


Other 




Ref. 


Spectral 








kms -1 




Type 


11 


68273 


7 2 Vel 


1550 


c 


WC8+07.5III 


90 


156385 




2045 


a 


WC7 


135 


192103 


V1042 Cyg 


1405 


a 


WC8 


116 




HM19-3 


2700 


b 


WC5+08 



Table 2. Journal of optical and IR spectroscopic observations. 
All SWS observations were obtained with the AOT06 scan. 



WR Epoch Telescope Wavelength Exposure 

Inst. Range (fJ.va) Time (s) 



11 17 May 1996 ISO-SWS 2.38-45.2 7,876 

90 9 Mar 1998 AAT-RGO 0.50-1.03 5 

16 Feb 1997 ISO-SWS 2.60-19.6 10,068 

135 Sep 1991 INT-IDS 0.38-0.73 24 

20 Aug 1994 UKIRT-CGS4 1.03-1.13 64 
19-21 Aug 1994 UKIRT-CGS4 1.61-2.21 576 

21 Aug 1994 UKIRT-CGS4 2.30-2.51 128 

11 Nov 1996 ISO-SWS 2.38-45.2 6,538 
146 21 Jul 1996 INT-IDS 0.36-0.68 800 

19 Aug 1994 UKIRT-CGS4 1.03-1.13 64 

20 Aug 1994 UKIRT-CGS4 1.23-1.33 64 
19-21 Aug 1994 UKIRT-CGS4 1.61-2.21 224 

21 Aug 1994 UKIRT-CGS4 2.30-2.51 64 

12 May 1996 ISO-SWS 2.60-19.6 16,922 



et al. 1990; St Louis et al. 1993), except for WR146 for which 
a terminal velocity was obtained in Paper I from ISO/SWS 
spectroscopy. In this section, we discuss the ground-based 
(AAT, INT, ESO) and space-based (ISO, IUE) observations, 
the journal of which is presented in Table |j| 

2.1 Ultraviolet spectroscopy 

Our UV dataset was obtained solely from the Interna- 
tional Ultraviolet Explorer (IUE) archive. Available high 
dispersion (HIRES), large aperture, short-(SWP) and long- 
wavelength (LWP, LWR) observations of WR90 and WR135 
were combined to provide high S/N datasets (St Louis 1990). 
Absolute flux calibration was achieved following the calibra- 
tion curve obtained by Howarth & Phillips (1986). Agree- 
ment between our final HIRES datasets and archival low 
resolution (LORES) observations was found to be excellent. 

Willis et al. (1986) have discussed the UV spectral mor- 
phology of WC stars, including WR90 and WR135. The 
principal spectral features for both stars are Siiv AA1393- 
1407, Civ AA1548-51, Hen A1640, and Cm A2297. A forest 
of Fev-vi lines are present in both stars, with FelV also 
prominent in WR135. 

2.2 Optical spectroscopy 

Previously unpublished spectroscopic observations of WR90 
and WR135 were obtained at the 2.5m Isaac Newton Tele- 
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scope (INT) and 3.9m Anglo Australian Telescope (AAT). 
The INT dataset was obtained with the IDS, 500mm cam- 
era, and GEC CCD during 1991 September, while the AAT 
dataset was obtained with the RGO spectrograph, 25cm 
camera, and a MIT/LL CCD during 1998 March. Each 
dataset was obtained using CCDs and reduced in a stan- 
dard manner using software available on STARLINK. Spec- 
tra were debiased, flat-fielded, optimally extracted, wave- 
length calibrated, and flux calibrated in the case of the INT 
dataset. Flux calibration was achieved for the AAT observa- 
tions via scaling to the level of the blue CTIO spectropho- 
tometry taken from Torres-Dodgen & Massey (1988). 

Arc spectra provided a measure of the instrumental res- 
olution: 2-3A (INT) and 4A (AAT). Complete spectral cov- 
erage between 3800-7300A (INT) and 5000-10300A (AAT) 
was achieved with, respectively, six and one grating settings. 
The atmospheric absorption bands removed using suitable 
comparison stars. 

The optical spectral appearance of WR146 has pre- 
viously been discussed in Paper I, while Dougherty et al. 
(2000) present a new high quality blue optical observation 
of that star. The morphologies of WR90 and WR135 are 
relatively similar, dominated by He 1-11 and C ni-iv emission 
features, except that the emission lines of WR90 are sub- 
stantially broader. C ill AA4747-51 and C IV AA5801-12 are 
the two strongest features in each case. Several O ni-v fea- 
tures are present around AA2950-3150, and A5592 with Cll 
present at A4267 for WR135. 

2.3 Near-IR spectroscopy 

Our principal near-IR dataset was obtained at the 3.8m 
U.K. Infrared Telescope (UKIRT) with the cooled grating 
spectrograph CGS4, the 300mm camera, a 751/mm grat- 
ing and a 62x58 InSb array in 1994 August. Observations 
of WR135 and WR146 were bias-corrected, flat-fielded, ex- 
tracted and sky-subtracted using CGS4dr (Daly & Beard 
1992). Subsequent reductions and analysis were carried out 
using FIGARO (Shortridge et al. 1999) and DIPSO (Howarth 
et al. 1998). In order to remove atmospheric features, the 
observations were divided by an appropriate standard star 
(whose spectral features were artificially removed) observed 
at around the same time and similar air mass. In regions 
of low atmospheric transmission at UKIRT the reliability 
of line shape and strength must be treated with caution 
(e.g. the Pa region) . The near-IR morphology of WR135 and 
WR146 have been discussed previously by Eenens, Williams 
& Wade (1991). 

We also utilise intermediate dispersion CCD spectra 
of WR135 covering 0.97-1.03/im, published by Howarth & 
Schmutz (1992), and obtained at the INT with the Inter- 
mediate Dispersion Spectrograph (IDS) and a GEC CCD 
in 1990 October. Although Howarth & Schmutz (1992) ob- 
tained observations of the He I AIO83O+P7 line, we prefer 
to use our new, lower resolution UKIRT observations of this 
feature. (The sharp fall in efficiency of the GEC CCD long- 
ward of P7 leads to an ill-defined continuum.). We attempt 
an approximate flux calibration for these data, by setting 
the local continuum via interpolation of the observed 0.7/im 
(INT) and 1.03^m (UKIRT) flux levels. 



Table 3. Summary of photometry, reddenings and distances for 
programme stars, including comparisons with reddenings from the 
literature (lit). Observed magnitudes correspond to the systemic 
value for binary systems, and include the Schmutz & Vacca (1991) 
corrections to Smith (1968) photometry (*) 



WR 


„sys 


Ref. 


E B -v (lit) 


R 


Ref. 


M„ wc 


d 


Ref 




mag 




mag 






mag 


kpc 




11 


1.70* 


a 


0.04 (0.03) 


3.1 


e 


-3.7 


0.26 


h 


90 


7.41* 


a 


0.38 (0.44) 


3.1 


d 


-4.7 


1.55 


e 


135 


8.51 


b 


0.37 (0.35) 


3.1 


a 


-4.3 


2.09 


f 


116 


13.91 


c 


2.87 (2.80) 


2.9 


c 


-5.3 


1.40 


g 



a: Smith (1968), b: Massey (1984), c: Paper I, d: Morris et al. 
(1993) e: Schaerer et al. (1998), f: Lundstrom & Stenholm 
(1984), g: Dougherty et al. (2000), h: van der Hucht et al. (1997) 



2.4 Mid-IR spectroscopy 

New mid-IR data of WR90 and WR146 were obtained 
as part of Guest Observer programme, AJWWOLF (P.I. 
Willis), with the Short Wavelength Spectrograph (SWS; de 
Graauw et al. 1996) onboard the ESA Infrared Space Ob- 
servatory (Kessler et al. 1996). Observations of HD 117297 
(WR53), HD 164270 (WR103) and HD 165763 (WR111) 
from the AJWWOLF programme were of insufficient qual- 
ity to provide mid-IR line flux measurements, and so these 
were excluded from the present study. Consequently, we 
have included SWS observations of WR11 and WR135 from 
the Guaranteed Time programme WRSTARS (P.I. van der 
Hucht). In all cases the SWS AOT6 observing mode was 
used to achieve full grating resolution, A/ A A ~ 1300 — 
2500. The continuous wavelength coverage was 2.60-19.6 /im 
for datasets from the AJWWOLF programme, and 2.38- 
45.0^m for WR135 and WR11. 

Data of detector "bands" 1 and 2 respectively cover 
wavelengths of 2.38-4.08 /an using 12 In:Sb detectors and 
4.00-12.05 /im with 12 Si:Ga detectors, employing entrance 
slits that give an effective aperture area of 14 x 20 arcsec 
on the sky. Band 3A to 3D data cover 12.0-27.6 /im using 
12 Si:As detectors, with sky coverage of 14 x 27 arcsec. fi- 
nally, band 3E and 4 data cover 27.5-45. O^rn using 12 Ge:Be 
detectors, with sky coverage of 20x33 arcsec. 

The total integration time was set to allow one complete 
scan over the wavelengths selected within each "AOT band" , 
defined by the permissible combinations of detector band, 
aperture, and spectral order (cf. de Graauw et al. 1996). The 
observation time includes dark current measurements, and a 
monitor of photometric drift for the detectors of bands 2-3. 
A drift measurement is not normally made for the relatively 
stable In:Sb detectors. 

The data processing of the SWS data for WR146 was 
discussed in Paper I, with WR11, WR90 and WR135 re- 
duced in a similar manner. The stellar spectrum of WR11 
will be discussed elsewhere (Morris et al. 2000, in prep), with 
solely mid-IR fine structure line fluxes measured here. 
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3 INTERSTELLAR REDDENING AND 
DISTANCES 

In this section, interstellar reddenings and distances to the 
programme stars are discussed. Pre-empting results from 
Sect. reddenings are directly obtained from comparing 
theoretical synthetic spectra with de-reddened observations. 
Distances either follow from cluster/association member- 
ship, or assumed absolute magnitudes. Table ^ provides a 
summary of the derived reddenings and distances for our 
sample of stars. A comparison with reddenings from the 
recent literature shows good agreement. The distances to 
WR135 (Cyg OB3 member) and WR11 (hipparcos) are 
known with confidence, while those adopted here for WR90 
and WR146 deserve comment. 



3.1 WR90 

WR90 is not a member of an association or cluster. We 
therefore estimated its distance based on the mean abso- 
lute visual magnitude of other Galactic WC7 stars. Un- 
fortunately, all five WC7 stars that are members of as- 
sociations/clusters are within binary systems (Lundstrom 
& Stenholm 1984). We restricted the sample to those for 
which reliable reddenings were known (HD97152=WR42, 
HD152270=WR79 and HD192641=WR 137) from Morris et 
al. (1993). An absolute visual magnitude of M v — —4.7 
mag(cr=0.8) was obtained for the WC7 components by com- 
paring their emission line strengths with WR90 (specifically 
Ciii-iv A4650, CivA5804). Our estimate is in good agree- 
ment with van der Hucht et al. (1988) and Smith et al. 
(1990) who derived M v — —4.8 mag. Pre-empting results 
from Sect. ^ we derive Eb-v=0-38 mag, so that the uncer- 
tainty in the calibration and reddening (± 0.02 mag) implies 
d=1.55kpct j[. 



3.2 WR 146 

Dougherty et al. (1996) used IR and mm photometry to 
estimate a distance of 1.2±0.3kpc towards WR146. Radio 
observations revealed two separate components, namely the 
(thermal) WC emission plus the (non-thermal) bow shock 
emission between the two winds. 

Meanwhile, a lower distance of 0.75±0.15kpc was de- 
rived in Paper I, using the mean M v for a WCE star (Smith 
et al. 1990), plus an assumed spectral type of 08.5 V for 
the companion. Subsequently, Niemela et al. (1998) used 
WFPC2 aboard the Hubble Space Telescope (HST) to mea- 
sure the individual UBV magnitudes of the WC and OB 
components of WR146. They derived V=13.64 mag for the 
WC star, 0.24 mag brighter than the OB component. In 
contrast, a difference of 0.8 mag was assumed in Paper I. 

The combined HST and radio datasets indicate that 
the OB companion possesses a powerful stellar wind, with 
a giant or supergiant luminosity class. Recent optical spec- 
troscopy supports a supergiant 08 classification (Dougherty 
et al. 2000). As discussed by Dougherty et al., a large dis- 
tance to this system would result if the absolute magnitude 
of the companion was typical of late O supergiants (Conti & 
Alschuler 1971). Instead, we adopt the distance of 1.4kpc, as 
derived by Dougherty et al. At this distance, WR146 would 



Table 4. Summary of the WC model atoms, and source of atomic 
datasets. N F is the number of full levels, N s the number of super 
levels and N TRA ns the corresponding number of transitions. The 
last column refers to the upper level of a given ion included in 
our treatment 



Species 


iV F 


N s 


-^Vtrans 


Ref 


Details 


He I 


39 


27 


315 


a 


n <14 


He II 


30 


13 


135 


b 


n <30 


He III 


1 


1 








C II 


88 


39 


791 


c,d 


nl <2p3d 4 D° 


C in 


243 


99 


5513 


e,f 


nl <2p4d 1 D° 


C iv 


64 


49 


1446 


g 


n <30 


C v 


1 


1 








O ii 


3 


3 


3 


c.h 


nl <2p 3 2 P° 


O in 


50 


50 


213 


c.i 


nl <2p4f 1 D 


O iv 


72 


30 


835 


e,j 


nl <2p3p" 2 P° 


O v 


91 


31 


718 


e,k 


nl <2p4p 3 P 


O vi 


19 


13 


72 


g 


n <5 


O vn 


1 


1 








Si iv 


28 


17 


129 


h 


n <6 


Si v 


1 


1 








Fe iv 


280 


21 


5055 


1 


nl <3d 4 ( 1 G)4p 2 P° 


Fe v 


182 


19 


2517 


m 


nl <3d 3 ( 2 D)4p 1 P° 


Fe vi 


80 


10 


722 


n 


nl <3d 2 ( 1 S)4p 2 P° 


Fe vn 


153 


11 


1213 


o 


nl <3p 5 ( 2 P)3d3(b 2 D 1 P° 


Fe VIII 


1 


1 










1299 


440 


20007 







(a) Fernley et al. (1987); (b) Wiese et al. (1966); (c) Yu Yan 
et al. (1987); (d) Yu Yan & Seaton (1987); (e) Nussbaumer & 
Storey (1983, 1984); (f) P.J. Storey (unpublished); (g) Peach et 
al. (1988); (h) Seaton (1995); (i) Luo et al. (1989); (j) Luo & Prad- 
han (1989); (k) Tully et al. (1990) (1) Becker k, Butler (1995b); 
(m) Becker & Butler (1992); (n) Becker & Butler (1995a); (o) 
K. Butler (unpublished) 

be a foreground object to Cyg OB2 (Lundstrom & Stenholm 
1984; Torres-Dodgen et al. 1991). 

We derived -Eb-v=2.87 mag and R=2.9 using the 
optical-IR reddening law of Steenman & The (1989, 1991). 
This provided a superior comparison between theoretical 
predictions and de-reddened observations than with either 
the Howarth (1983) or Cardelli et al. (1989) laws. Therefore, 
M v — —5.3 mag for the WCE component, by far in excess 
of 'typical' WCE stars (Smith et al. 1990). 



4 MODELLING TECHNIQUE 

Before discussing the fine analysis of each programme star, 
we introduce the spectroscopic technique followed here. We 
use the non-LTE code of Hillier & Miller (1998) which 
iteratively solves the transfer equation in the co-moving 
frame subject to statistical and radiative equilibria in an 
expanding, spherically symmetric and steady-state atmo- 
sphere. Relative to earlier versions of this code (Hillier 1987, 
1990), two major enhancements have been incorporated, of 
particular relevance to WC-type stars, namely (i) line blan- 
keting, (ii) clumping. Specific details of the techniques used 
are provided by Hillier & Miller (1998, 1999), with only a 
brief overview given here. 

As discussed by Hillier & Miller (1998, 1999), extremely 
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complex atomic models are necessary for the quantitative 
analysis of WC stars, a computationally demanding require- 
ment. Consequently, we make use of the technique of 'super- 
levels', in which several atomic levels of similar energies and 
properties are combined into a single one, a super level, with 
the populations of the super level calculated in the rate equa- 
tions. Populations of individual atomic levels are then calcu- 
lated by assuming that it has the same departure coefficient 
as the corresponding super level to which it belongs. In this 
way, extremely complex atoms of helium, carbon, oxygen 
and iron can be considered. In Table ^ a total of 1299 full 
(TVf) atomic levels are combined to produce just 440 su- 
per levels (N a ). The number of iron transitions is limited 
by gf < 10~ 4 , so that a total of 20,007 transitions are in- 
cluded. In addition to the specific models discussed here, we 
have also considered cases with expanded atomic datasets 
and ionization stages (e.g. Fein), for which the differences 
in emergent synthetic spectrum were negligible. 

Oscillator strengths, collision and photoionization cross- 
sections are taken from a wide variety of sources (Table |^). 
The OPACITY project (Seaton 1987, 1995) formed the ba- 
sis of most radiative rates, supplemented by calculations for 
CNO by Nussbaumer & Storey (1983, 1984) and Storey (un- 
published), and for iron by Becker & Butler (1992, 1995ab) 
and Butler (unpublished). 

The stellar radius (R*) is defined as the inner boundary 
of the model atmosphere and is located at Rosseland optical 
depth of 10 with the stellar temperature (T*) defined by 
the usual Stefan-Boltzmann relation. Similarly, T 2 /3 relates 
to the radius (R2/3) at which the Rosseland optical depth 
equals 2/3. 

Although the majority of spectral analyses of WR stars 
have adopted a standard /3=1 velocity law, there is both ob- 
servational and theoretical evidence for a more slowly accel- 
erating outflow (e.g. Schmutz 1997; Lepine & Moffat 1999). 
Consequently, we adopt a form for the velocity law (Eqn 8 
from Hillier & Miller 1999) such that two exponents are con- 
sidered with the result that acceleration is modest at small 
radii, but continues to large distances, i.e. 

_ ^0 + («oo - ^oxt - vo)(l - R+frf 1 + v cxt (l - R+frY 

1 + (VVcore)(exp((i?* ~ r) / h cS )) 

Here v ex t is an intermediate terminal velocity, v CO re the 
core velocity (typically a few kms -1 ), vo is the photospheric 
velocity (typically 100.0 kms -1 ), Voo is the final wind ve- 
locity and h c ff the scale height (~0.01ii*). For all models, 
we take /3i=l and /3 2 =50 as in Hillier & Miller (1999). 

There is now overwhelming evidence for the clumped 
nature of WR stars (e.g., Moffat et al. 1988; Moffat 1999), 
so we have adopted a simple filling factor approach. We as- 
sume that the wind is clumped with a volume filling factor, 
/, and that there is no inter clump material. Since radiation 
instabilities are not expected to be important in the inner 
wind we parameterise the filling factor so that it approaches 
unity at small velocities. Clumped and non-clumped spec- 
tra are very similar, except that line profiles are slightly 
narrower with weaker electron scattering wings in the for- 
mer. Although non-clumped models can be easily rejected, 
because of the severe line blending in WC winds M/y/J is 
derived by our spectroscopic analysis, rather than M and /. 
Also, this formulation introduces a revision of the BRA88 
analytical formula, which will be addressed in Sect 0. 



5 SPECTROSCOPIC ANALYSIS OF 
INDIVIDUAL STARS 

In this section, we discuss the analysis of our programme 
WC stars. Stellar parameters (T*, logL/L , M/y/f, C/He, 
O/He) were adjusted until the observed ionization balance, 
line strengths and de-reddened optical continuum flux dis- 
tribution were reproduced. Because of the substantial effect 
that differing mass-loss rates, temperatures and elemental 
abundances have on the emergent spectrum, this was an 
iterative process, in which initial model parameters were 
adopted from the literature if available (e.g. Koesterke & 
Hamann 1995; Morris et al. 1993). Terminal wind velocities 
are tabulated in Table |l| 

The wind ionization balance was ideally selected on 
the basis of isolated optical lines from adjacent ionization 
stages of carbon and/or helium. In practice this was diffi- 
cult to achieve because of the severe blending in WC winds. 
The low wind velocity of WR135 allowed He I A5876 and 
Hell A5412 to be used as ionization balance constraints. 
In other cases, the wind ionization balance was selected on 
the basis of carbon diagnostics. Unfortunately, from all the 
possible WC diagnostics, the usual classification lines C ill 
A5696/C IV A5804 are difficult to match, being so sensitive to 
minor stellar parameter changes. Consequently alternative 
diagnostics were sought. For WR90, our primary diagnos- 
tics were Cm A8500/Civ A7736 since these suffered from 
negligible contamination and are predicted by the model to 
vary smoothly across the temperature space. In the case of 
WR146 this spectral region, as well as lines of low ionization, 
were not available. We selected Hen A5412/Hei A10830, and 
Cm A6740/Civ 1.74/xm as our primary ionization diagnos- 
tics. 

Our experience with a variety of stellar models led us to 
select the strong ultraviolet P Cygni profiles of carbon (C ill 
A1909, 2297) as our principal mass-loss diagnostics. A major 
limitation with C IV A1550 was that different mass-loss rates 
also affected nearby Fe lines, which strongly modulated the 
predicted strength of the emergent P Cygni profile. Since 
ultraviolet observations of WR146 are unavailable, we relied 
on He 1 A10830 as the principal mass-loss diagnostic in that 
case. 

As in other recent spectroscopic studies of WC stars, 
Hell A5412/Civ A5471 were selected as the diagnostics for 
C/He determinations since the relative strength of these fea- 
tures are insensitive to differences of temperature or mass- 
loss rate. As discussed by Hillier & Miller (1998), this spec- 
tral region also contains a number of additional weak fea- 
tures, common to all stars. More problematic is that mis- 
leading C/He ratios would be obtained if (i) a limited num- 
ber of Civ atomic levels were included; (ii) homogeneous 
models were adopted in which electron scattering wings were 
incorrectly predicted. 

Test calculations were also performed using recombina- 
tion theory. Comparisons of C/He with our more sophisti- 
cated results for WCE stars was found to be reasonable - 
we obtain C/He=0.08 by number for WR146 in this study, 
in accord with the value obtained by Eenens & Williams 
(1992) based on IR recombination lines. WCL stars are more 
problematic since recombination coefficients for C ill are not 
available. In addition, the fraction of recombined helium is 
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Figure 1. Comparison between de-reddened (£?b— v=0.38, R=3.1) spectrophotometry of HD 156385 = WR90 (WC7) obtained 
from IUE, CTIO, AAT and ISO (solid lines), and line-blanketed, clumped model predictions (dotted lines). Stellar parameters are 
M=2.5xlO _5 M Q yr- 1 , T*=71 kK, log L/Lq = 5.5, C/He=0.25 and O/He=0.03 by number. Flux units arc in crgs/cm 2 /s/A 

© 1999 RAS, MNRAS 000, 



Neon abundances in WC stars 7 



not straightforward to assess in such studies, so that helium 
abundances could be underestimated. 

Oxygen abundances were more difficult to constrain, as 
already discussed by Hillier & Miller (1999), with the princi- 
pal diagnostic region spanning AA2900-3500. Since this spec- 
tral region was absent for WR146, we adopted C/0=4 by 
number, as predicted by current stellar evolutionary mod- 
els. Note that our neon abundance determination in Sect. [?] 
is relatively insensitive to the precise oxygen content. Re- 
garding silicon and iron, we adopt solar abundances since 
all distances are <2 kpc. 

We now proceed to discuss individual stars in detail. 



5.1 HD 156385 (WR90) 

The principal datasets for WR90 comprised HIRES IUE 
ultraviolet spectroscopy and AAT/RGO spectrograph op- 
tical observations. Secondary datasets were the blue CTIO 
spectroscopy from Torres-Dodgen & Massey (1988) and the 
2.6-5/im ISO spectroscopy (longer wavelength data were of 
insufficient S /N to be used as stellar diagnostics) . The com- 
bined, UV-optical-IR flux calibrated dataset for WR90 al- 
lowed a well constrained reddening of _Eb-v=0.38 mag, in 
agreement with the determination by Morris et al. (1993) 
(see Table ||). A distance of 1.55 kpc is implied from our 
assumed WC7 Af„-calibration. 

The difficulty in identifying the stellar continuum in the 
rich emission line spectrum of WC stars makes rectification 
imprecise, as emphasised by Hillier & Miller (1998, 1999). 
Fig. |l] demonstrates the excellent agreement between the line 
and continuum distribution of the model spectrum (dotted 
lines) and observations (solid lines), and includes the true 
theoretical continuum distribution (dashed lines). The num- 
ber of line features that are poorly reproduced is small, and 
includes Cm A5696 (too weak), A9710 (too strong), plus 
Civ A1548-51, Siiv A1393-1402 (both too weak because 
of Fe-absorption) and Ovi AA3811-34 (too weak). The use 
of absolute fluxes appears to give poor results in the line 
strengths around C ill A2297 since the de-reddened and the- 
oretical continuum do not exactly match. 

Our analysis reveals stellar parameters of T*=71,000K, 
log(L/L Q )=5.5 and M /^=8xl0~ 5 Moyr" 1 . Adopting a 
volume filling factor of /=0.1 indicates log M/(Mq yr _1 ) = 
—4.6 and a wind performance number of ~8. Use of C IV 
A5471/Hen A5412 revealed C/He=0.25±0.05 by number. 
An oxygen abundance of O/He=0.03±0.01 by number was 
obtained by matching Oiv AA2916-26, AA3063-72, and 
A3560-63. 



5.2 HD 192103 (WR 135) 

A substantial observational dataset was available for the 
WR135 analysis, particularly high quality UV (HIRES 
IUE), optical (INT), near-IR (INT, UKIRT) and mid- 
IR (ISO) spectrophotometry. An interstellar reddening of 
_Z5b-v=0.37 was obtained for WR135, in accord with Smith 
et al. (1990) (see Table |), but 0.18 mag lower than the 
more recent determination of Morris et al. (1993). Member- 
ship of Cyg OB3 implied an absolute visual magnitude of 
M v = -4.3 mag. 

Dereddened spectroscopy of WR135 (solid lines) is com- 
pared to our final synthetic model (dotted lines) in Figure til 
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Figure 2. Mo del predictions for the WC8 star WR135, analysed 
in Section (Top): Line formation regions for a selection of 
lines. The parameter £ is related to the observed flux emitted in 
the corresponding line, and defined as in Hillier (1987). The inte- 
gral of § over dlog(r/i?«) is proportional to the equivalent width. 
(Middle): Wind ionization stratification for carbon (solid) and 
helium (dotted). Note that helium does not recombine to its neu- 
tral state in the outer wind. (Bottom): Radial dependence of 
electron temperature (kK, dot-dash), density (cm -3 , solid), and 
wind velocity (kms -1 , dotted) 
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Figure 3. Comparison between the de-reddened (Eg_v=0.37 and R=3.1) spectrophotometry of HD 192103 = WR135 (WC8) obtained 
from IUE, INT, UKIRT and ISO (solid lines) and line-blanketed, clumped model predictions (dotted lines). Stellar parameters are 
M=1.3xlO" 5 M yr _1 , T+=63 kK, logL/L = 5.20, C/He=0.13 and O/He=0.03 by number. Flux units are in ergs/cm 2 /s/A. 
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Agreement is overall excellent, even for the forest of iron 
lines in the UV. The relatively low wind velocity of WR135 
permits a greater number of individual diagnostics to be 
selected, including Hell A4686, He I A5876, so that the de- 
rived properties can be treated with confidence. Indeed, C ill 
A5696 and C IV A5804 are fairly well matched in this case. 
The strong Cm spectral features at A2297 and A9710 are 
predicted to be 20-30% too strong. Other notable model de- 
ficiencies include underestimating the strength of C II emis- 
sion at A4267, A9900, and Ovi AA3811-34, with Om A3130 
too strong. As with WR90, the use of absolute fluxes, ap- 
pears to give poor results in the line strengths around He II 
A2530 since the de-reddened and theoretical continuum do 
not exactly match. This discrepancy simply reflects the in- 
adequacy of the adopted reddening law, rather than a fun- 
damental flaw with current models. 

We obtain the following stellar parameters: T*=63kK, 
log(£/L )=5.2 and M/v7=3.8x 1O" 5 M yr" 1 . Adopting a 
volume filling factor of /=0.1 indicates log M/(Mq yr -1 ) = 
—4.9 and a wind performance number of ~8. Use of C IV 
A5471/Hen A5412 revealed C/He=0.13±0.03 by number. 
Once again, the oxygen abundance is poorly constrained, 
although O/He=0.03±0.01 provides a reasonable match to 
O iv AA2916-26 and AA3063-72, with O III A3127 too strong. 

In Fig. ^, we present the wind structure of the fi- 
nal WR135 synthetic model, including line formation re- 
gions, ionization balance etc. Similar relations are shown 
for WR111 (WC5) in Hillier & Miller (1999). 

5.3 WR146 

Our principal observational datasets for WR146 are identi- 
cal to those used in Paper I, namely INT (optical), UKIRT 
(near-IR) and ISO (2.6-5/im because of the low S/N at 
longer wavelengths). The difference in our approach is to 
derive stellar and chemical properties solely from spectral 
synthesis, rather than recombination theory and indepen- 



dent modelling of the continuum. As discussed in Sect. 3.2 



we follow the distance estimate of 1.4kpc from Dougherty et 
al. (2000), in the light of new observations from Niemela et 
al. (1998). 

We have included the spectral energy distribution of a 
08 supergiant in our synthesis, using that for HD 151804 
(08If) derived by Crowther & Bohannan (1997). The IR 
free-free excess of this model is somewhat greater than 
equivalent temperature Kurucz (1991) models for which the 
lowest gravity available is log g=4. Because of the contami- 
nation from the late O supergiant, the line spectrum of the 
WC component of WR146 is relatively weak, as shown in 
Fig. The flux level of the O companion, the contribution 
of which declines with increasing wavelength, is illustrated 
as a dashed line in Fig. [|. Comparison between de-reddened 
observations and our synthetic spectrum is excellent, with 
the notable exception of C iv A5804. 

Our final stellar parameters for the WC5 star were, 
T„=57 kK, logL/L = 5.7, and Af/ v / 7=lxlO" 4 M yr- :L . 
A filling factor of /=0.1 reproduced the red wing of Cm- 
IV A4650-He II A4686. The high wind velocity of WR146 
meant that Civ A5471/Hell A5412 were blended. Never- 
theless, we were able to constrain the carbon content, de- 
riving C/He=0.08±0.02 by number, a factor of two times 
lower than that estimated in Paper I from recombination 



line analysis, but now in accord with Eenens & Williams 
(1992). Oxygen is extremely difficult to measure in WR146 
since the usual near-UV diagnostics are unavailable, so that 
O/He=0.02±0.01 is adopted. 

Finally, should the absolute magnitude of WR146 be 
more typical of other WCE stars, namely M v =— 3.7 mag 
(Smith et al. 1990), what would be the effect on the de- 
rived stellar parameters? In this case, our reddening towards 
WR146 would imply a distance of 660 pc, such that the com- 
panion star would be extremely faint M v =— 3.5 mag, typ- 
ical of an early B dwarf. The stellar properties of WR146 
would be unchanged, except that log L/Lq = 5.0, and 
Af/ v / 7=3.3xlO" 5 M yr- 1 . In this case it would be diffi- 
cult to reconcile these properties with the HST/radio ob- 
servations of WR146 (Dougherty et al. 1996; Niemela et al. 
1998). 



6 SUMMARY OF SPECTROSCOPIC RESULTS 

In this section we discuss the results of our quantitative 
analyses, and make comparisons with previous studies. 



6.1 Stellar parameters of WC stars 

In Table we provide a summary of the derived prop- 
erties of our programme WC stars, including results for 
WR11 and WR111 from De Marco et al. (2000) and Hillier 
& Miller (1999), obtained using identical techniques. De- 
spite our sample spanning WC5 to WC8, there is no ob- 
vious trend between spectral type and stellar temperature. 
Indeed, the WC5 component of WR146 is found to have the 
lowest stellar temperature, although this is certainly a most 
unusual WCE star! The subtle behaviour of C ill 5696A and 
C iv 5808A from our modelling indicates that, unfortunately, 
we are unable to use these as probes of wind ionization for 
WC5-8 stars (while we cannot obtain a simultaneous fit to 
both lines, a relatively small change in parameters can lead 
to a fit of either line) . 

Similarly, although there is considerable spread in car- 
bon abundances, WR146 exhibits the lowest carbon mass 
fraction, in contrast to the predicted increase in C/He ratio 
at earlier spectral type. Koesterke & Hamann (1995) de- 
rived a broad range of C/He ratios at each spectral type for 
WC5-8 stars. 

How do the present line blanketed, clumped results 
compare with previous studies? Unfortunately, WR135 is 
the sole programme star that has been the subject of quanti- 
tative studies in the past. Eenens & Williams (1992) derived 
elemental abundances from IR recombination lines, estimat- 
ing C/He=0.12 by number, in excellent agreement with our 
determination of 0.13. Howarth & Schmutz (1992) used a 
pure helium non-LTE model analysis to investigate WR135 
based solely on 1/im spectroscopy. Since pure helium models 
are expected to be inadequate for WC analyses, Koesterke & 
Hamann (1995) considered both carbon and helium in their 
study of WC stars, including WR135. 

In Table |^ we compare results from our present analy- 
sis of WR135 with these previous studies, scaling their re- 
sults to our absolute visual magnitude. Wind velocities and 
M/y/J are in relatively good agreement, while the effect of 
including carbon and blanketing has a dramatic effect on 
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Table 5. Derived stellar properties for programme WC stars, including results obtained by De Marco et al. (2000) for WR11 and Hillier 
& Miller (1999) for WR111 for comparison 
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Figure 4. Comparison between the de-reddened (i?B-v=2.87 mag, R=2.9) spectrophotometry of WR146 (WC5+OB) obtained from 
INT, UKIRT and ISO (solid lines) and line-blanketed, clumped model predictions (dotted lines). The contribution of the late O companion 
is taken into consideration, using model calculations by Crowther & Bohannan (1997) — its continuum distribution is indicated with 
dashed-lines. WC stellar parameters are M=3.3xlO _5 M yr" 1 , T*=57 kK, logL/L = 5.7, C/He=0.08 and O/He=0.02 by number. 
Flux units are in ergs/cm 2 /s/A. 



the stellar temperatures and luminosities, such that the lu- 
minosity derived here is a factor of two times higher than 
Howarth & Schmutz (1992), who adopted a stellar temper- 
ature of 35,000K. 

Stellar parameters are in much better agreement with 



Koesterke & Hamann (1995). They compared specific line 
strengths with model grids at fixed C/He ratio, and chose 
simple model atoms of Hel-n and Cii-iv. Although spec- 
tral comparisons between model predictions and observa- 
tions were not presented by Koesterke & Hamann (1995), 
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the fit quality was judged to be poor, with large discrepan- 
cies for the C ill A5696 and A6740 lines. Therefore, we have 
greater confidence in our results since detailed UV, opti- 
cal and IR synthetic spectrophotometry compare favourably 
with observations. 

From Table ^, blanketing and clumping conspire to re- 
vise the wind performance number, M /voo / (L/ c), from 30 
in the study of Koesterke & Hamann (1995) to just 5 for 
WR135! For our entire sample, performance numbers are 
<10, with previous studies indicating values of up to 100 
(Howarth & Schmutz 1992; Koesterke & Hamann 1995). 

6.2 Evolutionary status 

The wide range in stellar luminosity of our sample of WC 
stars, L/ Lq — 10 5 to 10 5 ' 7 , implies a considerable range in 
current stellar masses. Following the mass-luminosity re- 
lation for hydrogen-free WR stars of Schaerer & Maeder 
(1992), present masses of 9M (WR135), 13M (WR90) 
and 19Mq (WR146) are implied. In contrast, initial mass 
estimates are much more dependent on specific evolution- 
ary tracks. 

All models predict stars to pass through the WC phase 
at ages of between 2.7-4.5Myr. For initial 40 and 60Af© 
models, stellar masses during the early WC phase are pre- 
dicted to be ~14M0 and ~24M©, respectively. At higher ini- 
tial mass, the situation is extremely model dependent. For 
example, an evolutionary model with initial mass 120M Q 
will produce a 60Mq WC star using standard de Jager et 
al. (1988) mass-loss rates. Alternatively, a lOM© WC star 
will result based on (2x) enhanced mass- loss for the post 
main-sequence and WNL phases (Schaller et al. 1992), or 
only a 4M WC star for (2x) enhanced mass- loss during 
the entire stellar evolution (Meynet et al. 1994). Therefore, 
the present masses and ages of our programme WC stars are 
well constrained, but initial masses are not. 

In order to attempt estimates of initial masses, we com- 
pare our derived (C+0)/He ratios versus stellar luminosi- 
ties with evolutionary predictions in Fig. ^ Predictions for 
4O-12OM0 initial models assume (2x) enhanced mass-loss 
relative to de Jager et al. (1988) during the entire stellar 
evolution (Meynet et al. 1994). These models favour 60AfQ 
for WR146 and WR90 and 4OM or 85M© for WR135. Con- 
sequently, multiple initial mass estimates may result for in- 
dividual stars. Nevertheless, we currently favour 4OM0 for 
WR135 and 60M Q for WR146 and WR90. 

Specific comparisons between the stellar properties of 
WR90, WR135 and WR146 and Meynet et al. (1994) evolu- 
tionary predictions are made in Table (ij. Both the initial 
6OM0 and 85Mq models adopt mass-loss rates that ex- 
ceed observations by large factors during early WC stages. 
This is because evolutionary models for hydrogen-free WR 
stars adopt mass-loss rates that are solely functions of mass 
(Langer 1989). Clearly, future evolutionary models should 
take allowance for more appropriate WR mass-loss rates. 

6.3 Radio fluxes 

We now compare predicted radio fluxes from our spectro- 
scopic mass-loss rates. Hillier & Miller (1999) provide a for- 
mulation based on Eqn 9 in Wright & Barlow (1975), allow- 
ing for the filling factor according to Abbott et al. (1981). 
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Figure 5. Comparison of measured (C+0)/He ratios versus lu- 
minosities for our programme stars, including WR111 (WC5) 
from Hillier & Miller (1999), with theoretical expectations from 
the solar metallicity evolutionary tracks assuming enhanced (2x) 
mass-loss (Meynet et al. 1994) 



Table 6. Comparison of our derived stellar parameters for 
WR135 with Howarth & Schmutz (1992, HS92) and Koesterke 
& Hamann (1995, KH95). We have adjusted their parameters to 
our assumed visual absolute magnitude (M v = —4.3 mag) 
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In all cases, the dominant ionization at the outer bound- 
ary of our models, N e ~ 10 s cm~ 3 or 200i?* is He+, C 2+ and 
2+ (see Fig.g). Since the radio emitting region lies at lower 
densities, we also consider cases in which C + and/or + are 
dominant. Outer wind electron temperatures are typically 
8000K (Fig. |), except for WR146 for which the lower cool- 
ing produced by less metals implies 9000K. 

We compare predicted radio fluxes with observed val- 
ues taken from the literature in Table ^, in which a uniform 
UV /optical filling factor of /=0.1 has been adopted through- 
out. Note that the quoted radio flux for WR146 refers solely 
to the WC component (the southern mean emitted flux S5 
in Dougherty et al. 2000) . We find that consistency is excel- 
lent for WR135 and WR146 for doubly ionized carbon and 
oxygen, while singly ionized carbon and oxygen are favoured 
for WR90. However, C + is not predicted to be the dominant 
ionization stage in the outer wind of these stars (Fig. ^) , so 
that the predicted radio flux for this star appears to be too 
high. 

It is possible that the filling factor in the radio and 
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Table 7. Comparison between our derived stellar parameters for 
WC stars and evolutionary predictions from Meynet et al. (1994) 
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Table 8. Comparison between predicted and observed 6cm 
(4.9GHz) radio fluxes in our programme stars, assuming /=0.1 
for the UV/optical line forming region. The filling factor entry in 
the table corresponds to the adopted value in the radio region, 
indicating the change in that quantity required to recover the ob- 
served flux at 4.9GHz. We have included WR11 in this table, since 
a neon abun dance determination for this star is to be carried out 



in Sect. 7.2 
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optical forming regions differs significantly. The predicted 
radio flux of WR90 would be in excellent agreement with 
the observed value if the radio filling factor was ~0.2, while 
that of WR11 requires ~0.05. Hillier & Miller (1999) found 
a similar discrepancy in their study of WR111. Note also 
that there is observational evidence that WR90 may not be 
single, since it is a non-thermal emitter (Leitherer et al 1997; 
Chapman et al 1999). 

We also include calculations for WR11 in Table ^, since 
we attempt to re-derive its neon abundance determination 



in Section (?]. As for the other stars in our sample, possible 
variations in filling factor between the inner and outer wind 
are important. 



7 NEON AND SULPHUR ABUNDANCES IN 
WC STARS 

We are now in a position to determine neon abundances 
in our programme WC stars based on ISO spectroscopy, 
supplemented by sulphur determinations for WRU. In this 
section we first provide a revised formulation for the deter- 
mination of neon in a clumped medium, following BRA88, 
and subsequently provide measurements for each star. 



7.1 Ionic abundances from fine-structure lines in 
an inhomogeneous wind 

We have re-derived the numerical and analytical forms for 
the determination of ionic abundances in clumped winds 
from BRA88. We consider a fine-structure line from ion i, 
with transition energy hv u i. If D is the distance to the star 
and I u i is the observed line flux, then 



n u A u ihv u iknr fdr 



erg s 



(1) 



where A u i is the line transition probability and we have in- 
troduced the filling factor / into their formulation. n u rep- 
resents the density of ions in the upper level, and can be 
written as, 



flu — fu^i 



(2) 



where rn is the species ion density and f u is the fractional 
population of the upper level. Upper level populations, /„, 
were determined for each ion by solving the equations of 
statistical equilibrium using equib (Adams & Howarth, priv. 
comm) for >30 electron densities in the range 10° to 10 12 
cm -3 , and 9 electron temperatures covering 5k to 14kK. 
Thus, 



fr 2 



(3) 



where 74 is the fraction of all ions represented by ion species 



7i 



(4) 



and A is the mass loss parameter (Eqn 8 from BRA88). 
Combining equations |^-^, the filling factor term cancels out, 
leaving 



7i 

Iui = -rr^A u ihv u iA 



f u (r,f,T)dr 



erg cm s 



"(5) 



We deviate from BRA88 by carrying out the integral in den- 
sity, rather than radial space. Since r(N e ) is a bijection, we 
can modify this integral so that the dependency of f u on 
A and / (and hence M) is removed from the integral term. 
Equation H can then be modified to: 



Ur.1.T),lr= J^j- 



MN e ,T) 

m- 6 



(6) 
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Figure 6. ISO/SWS observations of fine structure [Siv] 10.5^im, 
[Ncn] 12.8/im and [Ne III] 15.5/mi emission in WRII, WR90, 
WR135 and WR146. 



For a more straightforward numerical computation, Equa- 
tion ^| is slightly adjusted: 



MN e ,T) 



dN e = Zn(10) 



fn{N e ,T) 



dlog{N e ) (7) 



The effect of this modification leads to a significant improve- 
ment over BRA88, who used an average value for each re- 
gion Ar. since the line formation region is a very sensitive 
function of { u . Also, the dependence of the integral on M is 
removed. Consequently, BRA88 overestimated the integral 
term, producing a lower elemental abundance. Our final nu- 
merical expression for the ion number fraction 74, is (cgs 
units) : 



7» = 



(47r^m ff u 00 ) 1,5 
Zn(10)/ ' 25 



M 



F U {T) y/^ylAuihvui 



with 

Fu(T) = 



fu(N e ,T) 
y/K 



dlog(N e ). 



(8) 



(9) 



Provided radio and recombination processes are used to de- 
rive M/y/J, elemental abundances are only weakly depen- 



dent on the distance (oc D 



since mass-loss rates depend 



on D . Therefore elemental abundances obtained with fill- 
ing factors of /=0.1 or 1.0 differ by a factor of 1.8. 

In order to allow for the possibility of a clumped 
medium, the analytical expression for the ion number frac- 
tion given by BRA88 (their equation A13) also needs to be 
multiplied by a factor of \[J . We include determinations of 
7; using both the integral and analytical expressions in our 



subsequent calculations. We find that the analytic expres- 
sion is reliable to within ~20%. 



7.2 Neon abundances for WC stars 

We are now in a position to evaluate neon abundances in 
our programme WC stars, using Equation |§| Figure ^ shows 
ISO/SWS observations of [Neni] 15.55/im emission for each 
of our programme stars, together with [Nelll 12.8/^m in 
WRII. Emission line fluxes are listed in Table H. 

Table ^| also includes atomic quantities used to derive 
ion fractions. Transition probabilities are taken from the 
NIST Atomic Spectra Database, while collision strengths 
are obtained from papers in the IRON Project series (see 
Table |^ for references). ISO flux measurements for WRII 
compare closely with previous measurements from IRAS and 
ground-based observations (van der Hucht & Olnon (1985); 
BRA88), as indicated in Table §. 

In the following comparison, we include abundances de- 
rived from ISO mid-IR neon lines for WRII, using stellar 
parameters derived by De Marco et al. (2000). 

Neon abundance determinations are sensitive to the 
ionization balance in the line forming region for the fine- 
structure lines, and is comparable to that of the radio form- 
ing continuum, iV e ~10 5 cm -3 . Ne + , with an ionization po- 
tential that is higher than C + and + , is observed solely 
in WRII (WC8). Therefore, we have allowed for the pos- 
sibility that the carbon and oxygen ionization balance are 
singly ionized for WC8 stars, with He + and doubly ionized 
carbon and oxygen otherwise. Table |l^ shows predicted ion 
abundances of Ne 2+ for each case, plus measured Ne + for 
WRII, with upper limits otherwise. In all cases Ne/He~3- 
4xl0 -3 by number, significantly greater than the expected 
cosmic value of Ne/He~5xl0 -4 in the C and O-enriched 
WC environment. 

It is possible that neon exists in (unseen) higher ioniza- 
tion stages, specifically Ne 3+ . However, the ionization poten- 
tial for Ne 3+ , 97eV, is significantly higher than C 3+ (64eV) 
and 3+ (77eV) which are not expected to be present in the 
outer winds of WC stars (Fig. |^). 

In summary, measurement of fine-structure lines of 
neon from ISO/SWS observations reveal Ne/He=0.003- 
0.004 (Ne/C~0.02), a factor of 6-8 times higher than cosmic 
abundances of Ne/He=0.0005 for the carbon rich winds of 
WC stars, supported also by sulphur abundance determina- 
tions. However, what additional sources of uncertainty are 
there in our Ne/He determinations, and why do our results 
for WRII differ from BRA88, who obtained Ne/He=0.001? 

Although the ISO-SWS neon line fluxes are in good 
agreement with the ground-based and IRAS neon line fluxes 
from BRA88 (Morris et al. 1998), we derive an elemental 
abundance that is three times higher. The source of this 
discrepancy is due to a different distance to WRII, and the 
use of a non-clumped mass-loss rate by BRA88. Morris et al. 
(2000) have recently emphasised the need for a reliable mass- 
loss rate estimate in the determination of neon abundances. 
In addition, our derived neon abundance is in good agree- 
ment with Morris et al. (1998) who combined ISO/SWS 
neon fluxes with the HIPPARCOS distance to WRII, and the 
X-ray derived mass-loss rate of Stevens et al. (1996). 

The other major factor affecting neon abundances is 
clumping. Our UV/optical analyses use filling factors of 
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f=0.1, yet we cannot observationally constrain the filling 
factor to better than 0.05< / <0.25, indicating a further 
20% uncertainty in Ne/He. Of greater importance, we as- 
sume identical filling factors for the optical line forming 
region (~10 n cm~ 3 ) and the neon emitting region (~10 5 
cm -3 ). Neon abundances would be increased by 40% for 
WR90, decreased by 30% for WR11, with WR146 and 
WR135 unchanged, assuming that (i) UV/optical mas s- los s 
rates are fully consistent with radio fluxes (recall Sect. 6.3), 



Table 10. Neon and sulphur abundances derived for the pro- 
gramme stars, using the stellar parameters given in the previous 
table. In all cases athe ionization balance is assumed to consist 
of He + , C 2+ and 2 + . For each star, the first entry refers to the 
7i derived from the (more reliable) integration method, with the 
second obtained from the analytical expression. 



and that (ii) volume filling factors in the neon emitting re- 
gion are identical to those in the radio region. 



7.3 Sulphur abundances for WR11 

In order to assess the reliability of our derived abundances, 
we have also calculated the sulphur abundance for WR11 
based on ISO observations of fine-structure [Siv] 10.5/im 
and [S in] 18.7/^m lines (Table ^|). Since sulphur is not en- 
hanced by nucleosynthesis, abundances should correspond 
to the cosmic value. The line forming region for the sulphur 
fine-structure lines peaks at iVe~10 4 cm -3 , somewhat lower 
than the neon lines. 

ISO spectroscopy of [S iv] 10.5/xm confirms the line flux 
measured from ground-based spectroscopy by BRA88, while 
stellar modelling anticipates weak contamination from lines 
of He i (12-8) 10.52/Ltm and Cm (20-17) 10.54/im. The stel- 
lar analysis of WR11 by De Marco et al. (2000) predicts 
a 15% contribution from these lines to the observed flux, 
which has been corrected accordingly. For the first time in a 
Wolf-Rayet star, ISO reveals the presence of [S ill] 18.7/im, 
blended with the stellar He I (14-10) 18.62/im feature, re- 
sulting in a 30% decrease in [S ill] flux. Use of an inappro- 
priate mass- loss rate and distance for WR11 led BRA88 to 
suggest S 3+ /He=2.5x 10~ 5 , such that they were obliged to 
predict S 2+ >S 3+ , with an expected high [S ill] 18.7/^m line 
flux that is not confirmed by ISO observations. 

From Table [ji], we find 7 S 2 + = 1.9 x 10~ 5 and 7 s3 + = 
3.1 x 10~ 5 , which imply S/He=6xl0~ 5 by number. This is 
in good agreement with the cosmic value of 7.5xlO _J for 
the C and He enriched environment of WR11. Therefore, 
our determinations imply Ne/S=50 for WR11, a factor of 
eight times greater than the cosmic value of Ne/S~7. Since 
these lines are formed in similar regions of the stellar wind, 
Ne/S abundances are essentially independent of clumping, 
and reveal a degree of neon enrichment relative to sulphur 
that is equivalent to that derived earlier for He. 



7.4 Comparison of abundances with theoretical 
predictions 

Regarding theoretical expectations, the level of neon enrich- 
ment is expected to be strongly correlated with the carbon 
content (Schaller et al. 1992; Meynet et al. 1994). In Fig.| we 
compare observed Ne/He and Ne/C versus C/He ratios with 
theoretical expectations from the Schaller et al. (1992) and 
Meynet et al. (1994) evolutionary models at solar metal- 
licity for an initial mass of 60Mq. We find that observed 
neon abundances are a factor of two below expectations. 
Specifically, we derive Ne/C=0.02±0.01 in all cases, while 
Ne/C>0.03 is predicted during this phase of the WC evolu- 
tion. Error bars shown in the figure account for uncertainties 
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in distance and (uniform) filling factors, but do not allow for 
the possibility of a varying volume filling factor between the 
neon and UV/optical line forming regions, which could be 
responsible for the discrepant cases. 



8 SUMMARY 

We have performed quantitative analyses of a small sample 
of WC5-8 stars, using models that account for line blanket- 
ing and clumping. Comparisons between synthetic spectra 
and de-reddened UV to mid-IR observations are excellent, 
with few modelling deficiencies identified. Stellar parameters 
support previous determinations (e.g. Koesterke & Hamann 
1995), except that the incorporation of blanketing yields 
higher stellar luminosities, while clumping indicates lower 
wind performance numbers, supporting the conclusions of 
Hillier & Miller (1999) for WR111 (WC5). Future studies 
will derive properties of WC-type stars, at both earlier (WO) 
and later (WC9) spectral type, and investigate whether pre- 
dicted ionizing properties are consistent with nebular obser- 
vations (e.g. Crowther et al. 1999). 

ISO/SWS spectroscopy reveals the presence of neon 
fine structure transitions, allowing abundance determina- 
tions. Using a revised formulation of the BRA88 technique 
to account for wind clumping, we derive neon abundances of 
Ne/He=3-4x 10 -3 by number, seven times higher than the 
cosmic value adjusted for the H-depleted WC environment, 
supported by Ne/S=50 for WR11 from sulphur fine struc- 
ture lines. The Ne enrichment is a factor of ~2 times lower 
than predictions of current theoretical models. However, dif- 
ferences in volume filling factors between the (high density) 
UV/ optical line formation regions and (low density) mid-IR 
fine-structure forming regions represent the greatest source 
of uncertainty in current Ne/He abundance determinations. 
Nevertheless, Ne/S provides an independent confirmation of 
the neon enrichment since it is independent of outer wind fill- 
ing factors. Future large ground and space-based telescopes 
that are optimised for the IR, will allow neon and sulphur 
line flux measurements and abundance determinations for 
more distant Wolf-Rayet stars. Of particularly interest are 
abundances in carbon and oxygen-rich WO stars. 
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Table 9. Observed mid-IR fine structure line intensities (units of 10 — 12 erg cm -2 s — 1 ), and adopted atomic parameters, including 
statistical weights of the upper and lower levels, u) u and oj;, transition probability, A u i, and collision strength Q u i at T e =8000K. Fluxes 
are obtained from ISO/SWS observations in all cases, except for literature measurements for WR11, obtained from IRAS/LRS (van der 
Hucht & Olnon 1985) and the UCL spectrometer at the AAT (BRA88). 
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(a) Naqvi (1951); (b) Galavis et al. (1995); (c) Saraph & Storey (1999); (d) Saraph & Tully (1994); (e) Butler & Zeippen (1994); 
(*) [Sm-iv] line fluxes are shown prior to correction for the presence of He I and Cm stellar features (see text). 
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Figure 7. Comparison between Ne/C (bottom) or Ne/He (top) 
ratios versus measured C/He, also confronted to theoretical 
expectations from the solar metallicity evolutionary tracks of 
Schaller et al. (1992) and Meynet et al. (1994) for 60M Q . Three 
tracks are shown: (i) normal mass-loss rates during the entire evo- 
lution of the star (dotted); (ii) enhanced (2x) mass-loss during 
the post main-sequence phase (dashed); (iii) enhanced (2x) mass- 
loss during the entire stellar evolution (solid). Error bars account 
for uncertainties in distance and filling factor (see text) 
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